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. The two-dimensional cage model for polymer motion is discussed with an emphasis on the effect of 

sideways motions, which cross the barriers imposed by the lattice. Using the Density Matrix Method 
as a solver of the Master Equation, the renewal time and the diffusion coefficient are calculated as a 
function of the strength of the barrier crossings. A strong crossover influence of the barrier crossings 
is found and it is analyzed in terms of effective exponents for a given chain length. The crossover 
04 ' scaling functions and the crossover scaling exponents are calculated. 
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(L) ■ I. INTRODUCTION 

The motion of a single polymer dissolved in a gel has been described in terms of reptation 0] . Typical for reptation 
d \ is that the polymer chain moves inside a tube, which can only be refreshed by growth and shrinkage at the ends. In 
^ ■ order to make this motion suitaWe for analysis, lattice models have been designed of which the cage model, proposed 
by Edwards and Evans in 1981 2], is the oldest. It was mainly considered as a simple model for reptation Q,y], which 
. is indeed the main mode of motion for polymers dissolved in a gel. The embedding lattice plays the role of the gel 
imposing barriers for the motion. The original introduction allowed for two types of motion: "reptations" , which are 
, motions along the confining tube and sideways motions or "barrier crossings" , in which the chain overcomes a barrier 
and thus changes the tube configuration. Sofar most attention has been paid to the reptations only. In this paper we 
concentrate on the interplay, which is quite delicate, as we will show. E.g. naively one might think that the diffusion 
coefficient is the linear sum of the contributions of the two mechanisms, but that is not at all the case, as has been 
noted by Klein Wolterink and Barkema Q in a similar context. This makes it difficult to analyze experiments, in 
which never the simultaneous presence of the two types of motion can be excluded. In the cage model the interplay 
of the two modes can be fully analyzed. 

The model has extensively been studied by Monte Carlo simulations, which seems the only way to deal with the 
QQ ' Master Equation for the stochastic motion [a, |^ 0, B 01 ■ Since the relevant Master Equation does not obey detailed 
PsJ . balance, no systematic solution method exists. The issues, to which we presently address ourselves, were the subject 
On of a lively debat in the early nineties. The various simulations showed certain tendencies, but the limitation to fairly 
short chains and the intrinsic statistical noise, prevented in our opinion definite settling the role of Rouse dynamics 
vs reptation. 

In this paper we use an alternative method, based on the analogy between the Master Equation and the Schrodinger 
, equation, by which the temporal evolution of the probability distribution of the chain configurations corresponds to 
■ the evolution of the wave function. Of course the wave function may be complex, while the probability distribution is 
^ , real and positive. Also the Master operator, viewed as a hamiltonian, is non-hermitian, which implies decay towards 
'"^ ■ the stationary state in contrast to the oscillatory temporal behavior of the eigenfunctions of quantum problems. 
^ Inspite of these differences one can benefit from the analogy, the more so because the Master operator corresponds 
O to the hamiltonian of a one-dimensional spin chain, for which the very efficient Density Matrix Method (DMRG) has 
been designed by White (IQ) . The cage model remains a one-dimensional quantum problem, irrespective the lattice 
in which it is embedded, because the chain itself is a linear structure. 

We focus on two dynamic properties: the renewal time r and the diffusion coefficient D and determine them directly 
from the Master Operator. Both properties refer to asymptotically long times (the stationary state) and thus our 
' calculations are complementary to the Monte Carlo simulations which probe the short and intermediate time behavior 
|lll | . The renewal time is the time needed for the chain to assume a new configuration, which has no memory of the 
original one. It is found from the gap in the spectrum of the Master Operator. The Master Equation always has 
a trivial eigenvalue 0, corresponding to the stationary state. Any other initial state ultimately decays towards the 
stationary state and the slowest relaxation time (the inverse of the gap) is the renewal time. Its calculation is difficult 
because for long chains, the gap gets very small and the excited states are hard to disentangle from the stationary 
state. In fact the gap decays with as a negative power z of the length N of the chain, such that r ^ . The zero 
field diffusion coefficient D is related to the drift velocity in a weak driving field. It decays as a power N^^ . 

We will confine ourselves to one- and two-dimensional embeddings. The degree of difficulty of the solution is related 
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to the dimension of the embedding lattice, simply because the higher the embedding dimension, the higher is the 
spin in the corresponding spin chain and the more states are required in the DMRG approximation. In this paper 
we will show that the one-dimensional version, which is admittedly unrealistic, allows an analytical solution. The 
two-dimensional embedding lattice presents the problem already in its full complexity, while of course the three- 
dimensional case is of the most experimental relevance. Actually the universal properties are believed to be the same 
for embedding lattices from d = 2 and higher. The practical limitation to two-dimensional embedding lattices derives 
from the fact that our computations are already at the limits of the present day possibilities. 
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FIG. 1: A picture of the polymer chain, consisting of reptons (grey dots) and some examples of allowed repton motions. The 
dashed arrows denote reptations (hernia moves), whereas the dotted ones represent barrier crossings. The solid arrow presents 
the driving field direction. 



The DMRG approach yields very accurate results in the domain where it converges. This enables us to use finite-size 
scaling analysis for the determination of the exponents and the crossover scaling functions. The case without crossing 
barriers (the reptation dynamics) has powers different from the case with these crossings (Rouse dynamics). From 
the viewpoint of dynamic scaling, the exponents are exotic and nice illustrations of how crossover takes places from 
one type of behavior to an other. A crossover scaling representation for r and D strongly elucidates their behavior. 



II. THE MODEL 



The model is a chain of + 1 reptons, connected by A^ links, (jji, • ■ • , i/n), to neighboring cells of a hypercubic 
lattice. A picture of a chain is given in Fig. ^ The mobile units are the reptons. The yi can take 2d values 
and because we have a two-dimensional embedding lattice, we can denote them by the directions A^(orth), £'(ast), 
5'(outh) and M^(est). Two consecutive reptons are always in adjacent cells, but the chain may backtrack and cells can 
be multiply occupied, (see Fig.^. The yi variables characterize the chain, since the absolute position is irrelevant for 
the properties that we consider. The statistics of the model is governed by the Master Equation for the probability 
distribution P{Y, t), where Y stands for the complete configuration (yi, • • • , yjv). It has the general form 

^^^^ = ^ [W^(Y|Y')P(Y', t) - W{Y'\Y)P{Y, t)] ^ ^ M(Y, Y')P{Y', t). (1) 



The Ws are the transitions rates of the possible motions that we have indicated in Fig. ^ The matrix M combines 
the gain terms (in the off-diagonal elements) and the loss terms (on the diagonal). M is a sum of matrices, for each 



repton one 



N 

M(Y,Y') = ^Af,(Y,Y'), (2) 

1=0 

where the sum runs over the reptons starting with the tail repton i = to the head repton i = N. The tail repton 
matrix is diagonal in all the link variables except the first 

N 

Mo(Y,Y')^mo{yi\y[)l[5y^,y,. (3) 

The tail repton produces exclusively reptations. Each move of the tail repton can be seen as a combination of a 
withdrawal towards the cell of the next repton and from thereon a move to a new cell. The matrix mg is explicitly 
given by the scheme 
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The parameter B = exp(e/2) is a bias, which accounts for the influence of a driving field, which can be an electric 
field when the reptons are charged. The value of (the small) e is a dimensionless measure for the strength of this 
driving field. The driving field is along the body diagonal, here in the North-East direction. So if the link N moves 
to the direction W , the tail repton moves two units in the direction of the field. The reverse process gets a bias B^^ . 
The head repton transition probabilities are given by a similar matrix with B^ replaced by B~^ . They are depicted 
in Fig. 121 One could give all the transitions an overall factor, but this would only influence the overall time rate. So 
we keep the unbiased transitions equal to 1. 

The internal repton i changes two consecutive links y^, yi+i- The matrix Mi is diagonal in all the other link variables. 
The transition matrix contains two types of contributions: 

• Reptations. These are the cases where y[ and y[j^i are opposite, e.g. N and S (sometimes called a hernia in the 
chain). Then repton i can retract to the cell to which it is doubly connected and then recreate a new hernia, 
a new pair of opposite links, e.g. E and W (see Fig.O. Note that the sequence EW differs from WE. Moves 
towards EW and WE get different biases. 

• Barrier crossings. A typical case is the sequence NE. It may flip to the sequence EN . To reach that new 
position it necessarily has to cross the lattice point enclosed by the two sequences. We give these transitions a 
(small) factor c together with the biases which measure the distance that repton i travels in the direction of the 
field (see Fig.lD). 

It is worth noticing that for the cage model all transition ratios are proportional to the bias B^ (or B~'^) [T^ . 

Generally, conservation of probability follows from the fact that the sum over the columns of the matrix M vanishes. 
So M has a zero eigenvalue and the eigenfunction corresponding to this eigenvalue is the stationary state of the 
system, to which every other initial state ultimately decays. The matrix is non-symmetric, due to the bias, which 
gives different rates to a process and its inverse. Thus one has to distinguish between left and right eigenfunctions. 
The left eigenfunction, belonging to the zero eigenvalue, is trivial (all components equal), the right eigenfunction is 
the stationary state distribution. 

The renewal time is usually defined for bias B — 1 (no driving field). For the diffusion coefficient we introduce a 
small driving field which induces an overall drift Vd ~ e. The proportionality coefficient gives the diffusion coefficient 
D according to the Einstein relation 
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FIG. 2: The allowed motions of the head repton with corresponding transition ratios. The dashed arrow presents the driving 
field. 



Af(Y, Y') is the matrix representation of the Master Operator Al, acting as a hamiltonian. As M is non-symmetric 
the hamiltonian M. is non-hermitian. One may view the states yi of the hnks as the states of a discrete plane rotator. 
However, translating the action of in terms of rotator operators, does not lead to a more transparent expression 
(except for d = 1, which maps on a spin 1/2 chain, see Section HV)l . 

The model has only a few parameters: the length of the chain iV, the strength of the driving field e and the relative 
strength c of the barrier crossings with respect to the reptating transitions. Experimentally the most interesting 
combination is the case where e — > and N ^ oo. This is a delicate limiting process since the product eN may stay 
finite and influences the nature of the stationary state strongly. The properties that we consider: renewal time r and 
diffusion coefficient D, refer to the case where this product remains infinitesimal. Thus effectively we have only N 
and c as parameters. We will see that also in this pair interesting scaling combinations occur. 



III. SYMMETRIES OF THE MASTER OPERATOR 



For our analysis optimal use of the symmetries of the Master Operator is vital. We have chosen the driving field 
in the North-East direction in order to make the directions N and E as well as S and W equivalent. We use two 
symmetry operations of the lattice: 

• Reflection with respect to the field axis. This turns the direction N into E and vice versa and similarly it 
interchanges S and W. We refer to this operation as S±. 

• Reflection parallel to the field axis, which interchanges the directions TV and W as well as E and 5". It is denoted 
as <S|| . 
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FIG. 3: The allowed motions of the hernia with corresponding transition ratios. The dashed arrow presents the driving field. 



The hamiltonian is invariant under the operation <S_l but not under It would, if it were accompanied by a field 
inversion. We can analyze the consequences of the symmetry by considering the following transformation 



(5) 



TZ is an orthogonal transformation of the states iV, E, S and W of the links. The new states have a definite symmetry 
character under the two operations. The state is even under S± and iSy , the state 1 is odd and even, 2 is even and 
odd and 3 is odd for both. Let us illustrate the effect of this rotation by applying it to the matrix mo, yielding 
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(6) 



The entries are now labeled by the states 0,1,2 and 3 and are given by 

mil = -2-(S2 + B-2) = m33, m22 = -2{B^ + B-^), 
m2o = -2(^2 _B-2) = 2m3i = 2mi3. 



(7) 



One sees that the matrix is block diagonal with a 2 x 2 matrix in the — 2 channel and one in 1 — 3 channel. This 
results from the invariance of the hamiltonian with respect to S± , since the states and 2 are even and the states 1 
and 3 are odd under this symmetry. So states of different parity under S± are not mixed by the hamiltonian. One 
also observes that the matrix becomes diagonal for B = 1. The off-diagonal elements concern transitions where the 
symmetry under field inversion is changed. 
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FIG. 4: The solid arrows present allowed motions for the barrier crossing with corresponding transition ratios. The dashed 
arrow presents the driving field. 



The symmetry character of the links in the new states can be carried over to larger segments of the chain, simply 
by multiplying the parities of the constituing links. The states 0,1,2 and 3 not only function as states of the links, 
but also as indices for the 4 different symmetry classes. What has been shown in detail for the tail repton, holds also 
for the total hamiltonian. It is invariant under and for i? = 1 also under iS|| . 

For the diffusion coefficient we expand the Master Equation in powers of e. 

M = Xo + eXi + ---, P(Y)=Po(Y)+ePi(Y) + --- (8) 

and obtain the equations 

XoPo = 0, A^oPi = -XiFo- (9) 

The first equation is trivially fulfilled by a constant Pq{Y\ since the matrix Mg is symmetric and the right eigenvector 
becomes equal to the trivial left eigenvector. Note that Po(Y) is the direct product of state for all the links. When 
we perform the rotation ||SJ| on all links, the vector Pq changes from a constant in all entries into the vector with a 1 
in the first entry and a in all others. 

The second equation is a set of homogeneous linear equations for the components of Pi(Y). It is soluble, since the 
right hand side of the equation is perpendicular to the left eigenvalue (which remains true for all orders in e) . So we 
can make the solution definite by requiring that it is also orthogonal to the trivial left eigenvector. Pi(Y) yields the 
lowest order drift velocity Vd- 

Prom the viewpoint of symmetries, the operator A4o is invariant with respect to both S\_ and . So it does not 
affect the symmetry character of Pi, which therefore inherits its symmetry from the right hand side of The latter 
derives its symmetry from M.\ since Pq is fully symmetric. From the example (jSJ we can see what the tail-repton 
hamiltonian does to the first link: it turns it from state (the building block of Pq) to state 2. Detailed calculations 
show that this also holds for the other components of A^i. By A^i the vector Pq turns from sector to sector 2. This 
is not surprising since closer inspection of the right hand side of shows that it is the microscopic expression for 
the drift velocity. This reverses sign under S\\ but stays invariant under 5^, which is indeed the symmetry character 
of sector 2. 
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FIG. 5: Comparison of m-dependence for the gap exponent at zero and non-zero c. 



IV. THE ONE-DIMENSIONAL MODEL 

As a curiosity we mention the cage model for a one-dimensional embedding lattice. Then the j/,; can take only two 
values, which one would naturally take as 1 and -1. 1 is a step forward and -1 a step backward along the embedding 
line. On the line there is no possibility of barrier crossing and we have therefore only reptations. An internal move 
takes place when a pair of links (1,-1) turns into (-1,1) or vice versa. In addition the tail and head link can change 
from 1 to -1 or from -1 to 1. Such a model belongs to the class of models which can be solved by the matrix product 
representation designed by Derrida et al. jlSj. Usually the model is discussed in terms of the variables 1 (a particle) 
and (a vacancy) and one visualizes the dynamics as a form of traffic. The particles cannot overtake each other and 
have to wait till they can exchange with a vacancy. The mutual exclusion of particles corresponds in the cage model 
to the fact that a repton cannot move when it is surrounded by two links of the same value. At the ends of the chain 
particles enter and leave with certain rates, which in the cage model means that the tail (and head) link can change 
into their opposite direction. Comparing the rules by which 1 and -1 can interchange in the cage model and the rates 
at which I's and -I's are created at the head and tail of the chain, one finds the rules for the equivalent traffic model. 

The class of such traffic models has been solved by Sasamoto 0| and independently by Blythc ct al. 15]. In 
this model has been related to the necklace model for reptation, where the results of the traffic model arc formulated in 
terms of the language of polymer motion and expressions are derived for the drift velocity and the diffusion coefficient. 
It is interesting that one has a solution for the whole range of values of the bias B. E.g. the value of the drift velocity 
becomes independent of the length N and reads for large N 

VdC:,^iB~B-'). (10) 

Thus indeed for small e the drift velocity becomes proportional to e. The diffusion coefficient follows using with 
an asymtotic decay as Z? ~ and the diffusion exponent equals a; = 1. 

Not only the gap can be calculated explicitly, but also the whole spectrum of the zero field haniiltonian, since it 
becomes equivalent with the Heisenberg ferromagnetic spin chain. The gap A reads 




A = -2(1 - cos(7r/7V)). 



(11) 
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FIG. 6: The comparison of m- dependence for the diffusion exponent at c = and c = 0.01. 



Note that the gap vanishes as A ^ N^"^ for long chains, yielding z — 2. 

Although all the moves of the d ~ 1 cage model are reptations, the exponents for the gap and the diffusion coefficient 
are not typical for reptation. The reason is that the one-dimensional model has no obstacles which slow down the 
motion by an order of magnitude in the chain length N, as is characteristic for reptation. 



V. EFFECTIVE EXPONENTS AND CONVERGENCE 



The DMRG expands the solution of the Master Equation in a basis of size m. The accuracy of the method can be 
tested by an internal parameter: the truncation error [T3 |. but more convincingly by the convergence of the results 
as function of m. A In-ln plot of the raw data for the gap as a function of chain length N is not very revealing, but 
rather misleading as has been pointed out by Carlon et al. |17| . A much more refined way of analyzing the data is to 
use effective exponents. For the gap we define 

_ lnr(Ar + 1) - lnT(7V - 1) _ dlnr 

" ln(iV + 1) - ln(iV - 1) ~ dhJ^' ' 

which is a function of the chain length. If the renewal time were a strict power law r ^ , the expression (|12|l would 
equal z, independent of N. We want to stress that only very smooth data can be used to calculate these effective 
exponents, because noise, which is e.g. inherent in simulations, will magnify in the ratio of small differences. 

In order to get an estimate of the convergence we give in Fig.Eltwo sets of curves for c = 0.01 and c — for various 
TO. The curves for c = 0.01 are quite close, such that further increase of m, does not lead to significantly different 
results. On the contrary, for the pure reptation case c = 0, the curves keep changing with to for longer and longer 
chains. Note that the basis is already unusually large (m = 180) for DMRG calculations. The large to was possible 
due to the speed-up of the process by using the full symmetry of the lattice. The double symmetry allowed to enlarge 
the size of the basis to by a factor 4. For each value of c there is a maximal value of A'' for which the result converge for 
feasible values of to. Within that range the DMRG values are also sufficiently accurate such that the small differences 
in H12() do not suffer from computational noise. In the pictures of the coming sections we only plot the data which do 
not depend on the value of m. 
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FIG. 7: The renewal time exponent as a function of the length of the chain for various values of the barrier crossing rate c. 
Note that we have shortened the curve for c = with respect to the best curve (m = 180) in Fig. |3 because for longer chains 
we have significant changes with increasing m. 



For the diffusion coefficient the domain of m-independent data is even more restricted. Here we introduce, similar 
to H12|l. the effective diffusion exponent 



xn 



\iiD{N +1) -liiD{N -1) ^ dlnD 
ln(Ar+ 1) -hi(iV- 1) ~ ^d\nN' 



(13) 



In Fig. 13 we show the exponent xjv for c = and c — 0.01. The latter is again reasonably convergent, but for the 
former we could not go to sufficiently large m such that a convergent domain starts to emerge. One also observes noise 
which is not visible on a In-ln plot, but which shows up as a result of small numbers in numerator and denominator 
in ifT^ . So practically our calculations are limited to values larger than c = 0.001. 

The local exponents zat for the renewal time obtained for various c are collected in Fig.[7| In Fig. |Slwe plot, in the 
same way, the local exponent xn for the diffusion coefficient. Generally, in spite of the above mentioned restrictions, 
we can make the following observations. 

• Chains of the order of ~ 100 are not yet in the asymptotic regime. The effective exponents still deviate 
appreciably from the asymptotic value. In other words, there are large corrections to scaling. In particular, the 
plateaus in the small c curves may easily lead to the conclusion that the exponent has settled on a too large 
value. 

• The influence of small values of c is quite strong for long chains. We come back on this point when we discuss 
the crossover behavior. 

• Although we have no clear evidence that the c = curve for the gap tends towards the asymptotic value Zqo = 3, 
it is clear that the curves for smaller and smaller c "try" to approach this theoretical value for reptation. The 
approach to the asymptotic value Zqo — 2, for larger c, is evident. This is the exponent for Rouse dynamics. 

• The curves that we could calculate for the diffusion coefflcient approach the asymptotic exponent Xoo — 1, which 
is again the Rouse exponent for diffusion. 
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FIG. 8; The diffusion exponent as a function of the length of the chain for various values of the barrier crossing rate c. 



VI. CROSSOVER SCALING 

As the curves of Figures and |H| show, the renewal time r and the diffusion coefRcient D are widely varying 
functions of the two parameters c and N. We can organize the data more transparently in terms of a crossover scaling 
function, aiming at data collapse. Anticipating the asymptotic values of the two regimes: c — > and a fixed c 7^ 0, 
the following representation is adequate for the renewal time 

t{N,c) ^ N^gic^N). (14) 

From such a representation one derives for the effective exponent the expression 

dlnr _ d\ng{c^N) 
dh[N ^ d\n{c^N) ■ 

The crossover function g{x) itself should be expandable for small arguments as 

9{x) = go+gix^ (16) 

and for large arguments as 

5(:^)=^-(5-i + — + ■••)• (17) 

a; V X J 

Inserting the asymptotic behavior (|17f) into (|14ll we obtain 

\^{T|N'^) = \^g^x-0\^^c^'■■■, (18) 

where the dots refer to corrections of order 1/A^. In Fig. |51we have made a plot of ln(r/A^^) vs Inc, extrapolated 
to A*" ^ c», which corresponds to the first two terms of H18|) . As one sees the curve is fairly straight, with a slope 
~ —0.51, in the domain where the data are most accurate. We use this value in a scaling plot of g(x), which is shown 
in Fig. 1101 The most important part of the figure is that for large argument the data collapse and fall on a curve 
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FIG. 9: Ln-ln plots of the renewal time and the diffusion coefficient as function of c. 



which decays as l/x, implying the crossover from reptation to Rouse dynamics. The pure reptation behavior foUows 
from the finiteness of g{x) for small argument. The value 5(0) can be derived from a plot of tN~^ versus N~^. We 
find the value g{0) ~ 0.026, which is in good agreement with the behavior of the curves for small argument. We see 
that the data at small x do not coincide as good as for large x. The reason is the inclusion of small values of N for 
small c. Here the N is not sufficiently large to have scaling. One should make N larger and c smaller to get good 
scaling in that region, but that regime is as yet inaccesible to us. 

In Fig. |51 also Inl/ DN is plotted as function of Inc. It gives again a straight line with the same slope 9 ~ 0.515. 
We use this value in the scaling plot. Fig. ^2 for the diffusion coefficient in the form 



D{N,c) = N-'^f{c^N). 



(19) 



As one sees the collapse is here amazingly good in view of the lesser quality of the data, as compared to the renewal 
time. Also smaller N and "large" c are included, although no real scaling can be expected for these values. It proves 
that crossover scaling works very well for the diffusion coefhcient. The crossover scaling function / approaches again 
a finite value at a; = that can be estimated from Fig. ^] as /(O) ~ 3.67. For large arguments, f{x) should behave 
as f{x —>■ 00) ~ X, which is confirmed by the plot. It shows again the crossover from reptation to Rouse dynamics. 



VII. DISCUSSION 



We have found that the renewal time and the diffusion coefficient can be transparently described by the crossover 
scaling functions H14|) and H19|) . In particular the diffusion data fit the scaling curve very well for practical all the 
points calculated. The crossover scaling exponent is found to be = 0.51 and we are fairly convinced that the exact 
value is 1/2. Not only the data support this value but also an analytical argument can be given for 9 = 1/2, which 
runs as follows. Reptation does not change the backbone of the chain (which results from stripping the hernias from 
the chain). The hernias walk along the backbone and are created and annihilated at the end of the chain. For the 
removal of a backbone segment of length N, the end repton has to diffuse over a distance of order N. For this diffusion 
along the backbone the curvilinear diffusion coefficient applies, which is an order N larger than the total diffusion 
coefficient, so it is of order (see also section llV)l . The time scale for diffusion is distance squared {N"^) divided 
by the diffusion constant {N~^). Thus a backbone segment of order N requires a time scale N'^/N~'^ = to be 
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FIG. 10: The crossover function g. 



renewed (this in fact explains the reptation exponent z = 3). On the other hand, to change that backbone segment by 
direct hopping over barriers, one needs the time N/c. The fastest process dominates and the competition is controlled 
by the ratio of the rates N^/{N/c) = cN"^. So the crossover scaling function should be a function of the ratio cN , 
which yields 9 — 1/2. 

One may wonder why it is so difficult for DMRG to reach long N for small c, while the hamiltonian simplifies 
for c — 0. As one observes from the figures [7| and |S1 there is a zone where the values make a turn from reptative 
behavior to Rouse dynamics. The basis in the DMRG approximation has to be large enough to notice this difference 
in behavior at the appropriate N, which grows as c~^/^. One needs a basis size m of the same order to describe the 
crossover. It does not mean that it is impossible to use DMRG for the pure cage model with c = 0. But then one 
has to use a representation of the hamiltonian which explicitly acknowledges the extra symmetries which are present 
in this model. Additional symmetries of the cage model are discussed in the papers of van Heukelum et al. 
similar situation occurs in crossover in the Rubinstein-Duke model 18]. 

The results of this paper are quite similar to the crossover found in a one- dimensional Rubinstein-Duke model with 
hernia creation and annihilation |l6j| . This leads us to believe that the crossover in gels is a universal phenomenon 
with the crossover exponent 9 = 1/2, independent of the embedding lattice, as long as the embedding lattice permits 
sideways motion, which cross the barriers. Sometimes the motion rules exclude crossing of barriers, e.g. in the one- 
dimensional embedding. But also in a lattice with triangular cells, crossing is impossible within the rule that links 
are always between nearest neighbor cells. 

We find that with fixed non-zero crossing rate c, the chain always tends towards Rouse dynamics for larger and 
larger N. This contrasts the general observation that in polymer melts the opposite tendency takes place: longer 
chains display reptative behavior ^| . It is clear that the obstruction due to other polymers cannot be seen as a 
fixed barrier, with a certain tunneling rate. Thus our results cannot be applied to polymer melts using a fixed rate for 
sideways motion. In other words, c must become a function of the chain length. To handle a chain length dependent 
hamiltonian gives a complication in DMRG. Apart from the difficulties to find an adequate model for polymer melts 
that allows to treat very long chains accurately, we may speculate that an "effective" rate c for sideways motion 
depends as a power N~" on the length N. The combination ciV^ ~ N'^~" determines whether one sees reptation 
of Rouse dynamics. It is tempting to take for a the effective renewal exponent z^, defined in (|12() . because it takes 
the renewal time for a polymer to get out of the way. As this a is always larger than 2, the combination shrinks 
with growing length, making the a even larger. Thus one observes the opposite crossover: from Rouse dynamics to 
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FIG. 11: The crossover function /. 



reptation. 
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